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Abstract 

Background: Quantitative PCR is rapidly becoming the standard method for analyzing gene expression in a wide 
variety of biological samples however it can suffer from significant error if stably expressed reference genes are 
not identified on which to base the analysis. Suitable reference genes for qPCR experiments on Staphylococcus 
pseudintermedius have yet to be identified. 

Results: Three reference genes in S. pseudintermedius were identified and validated from a set of eight potential 
genes iproC, gyrB, rpID, rho, rpoA, ftsZ, recA, sodA). Two strains of S. pseudintermedius were used, and primer 
specificity and efficiency were confirmed and measured. Ranking of the genes with respect to expression stability 
revealed gyrB, rho and recA as the best reference genes. This combination was used to quantify expression of a 
single biofilm associated gene, icaA, in logarithmic, stationary and biofilm growth phases, revealing that expression 
was significantly upregulated in the biofilm growth phase in both strains. 

Conclusion: Three reference genes, gyrB, rho and recA, were identified and validated for use as reference genes 
for quantitative PCR experiments in S. pseudintermedius. Also, the biofilm associated gene icaA was shown to be 
significantly upregulated in biofilm samples, consistent with its role in biofilm production. 
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Background 

Staphylococcus pseudintermedius is a common com- 
mensal organism of canines, but is also one of the most 
common causes of opportunistic infections [1-3]. Re- 
cently, methicillin resistant S. pseudintermedius (MRSP) 
has emerged and disseminated internationally [4,5], with 
2 major sequence types (ST68 and ST71) representing 
the variety of clinical infections in most regions [3]. One 
area of concern is the ability of this bacterium to produce 
biofilm something that might be an important virulence 
factor and complicate elimination of infections [6,7]. Ex- 
pression of genes pertaining to initial bacterial surface 
adherence and intercellular adhesion following biofilm 
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formation, such as microbial surface components recog- 
nizing adhesive matrix molecule (MSCRAMMs) which 
mediate cellular adhesion, and the intracellular adhesion 
{icaADBC) operon [8], reported to be at least partially re- 
sponsible for biofilm formation, likely affect the in vivo 
behavior of this organism, including resistance to therapy 
[9,10]. While there is significant postulation regarding 
these factors, understanding of the expression of anti- 
microbial resistance and biofilm associated genes in 
5. pseudintermedius and the subsequent clinical implica- 
tions is still poor. 

Quantitative real-time PCR (qPCR) is increasingly 
employed to quantify gene expression. While it can be very 
sensitive and specific, there are numerous pitfalls in its ap- 
plication that can easily result in misleading and incorrect 
conclusions. On of the most frequent errors is a failure to 
confirm the constitutive expression of the reference genes 
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used to measure the relative expression of genes of interest. 
Normalization of results in qPCR is vital to limit variability 
introduced by experimental conditions, sample preparation 
and analysis, and is one of the main underlying tenets of 
qPCR analysis. Selection of inappropriate reference genes 
can result in grossly incorrect conclusions owing to the 
miscalculation of gene expression. These and other reasons 
have prompted the development of minimum information 
for publication of qPCR experiments (MIQE) guidelines to 
ensure integrity, consistency and transparency of qPCR ex- 
periments, including standards for all aspects of experi- 
mental design, analysis and reporting [11]. 

Current recommendations suggest a minimum of 
three reference genes (ideally with M values below 1 
for heterogenous samples), and the inclusion of add- 
itional genes as necessary to obtain a pairwise variation 
value < 0.15 [12]. This requires specific validation of 
candidate reference genes in the bacterium (and ideally 
strains) to be studied, something that is lacking for 
S. pseudintermedius. While multiple qPCR studies have 
been performed in various Staphylococcus spp [13,14], 
study and validation of reference genes in S. pseudinter- 
medius is lacking, and it cannot be assumed that data 
from other staphylococci apply to this species. 

The objectives of this study were to evaluate several 
potential reference genes in S. pseudintermedius, to identify 
the optimum gene or gene combinations for future qPCR 
expression studies and to evaluate expression of icaA using 
validated reference genes. 

Results and discussion 

RNA isolation 

The modified protocol afforded good to excellent yields 
(mean yield 53.5 ± 23.4 ug, range 16.9 - 92.2 ug, concen- 
tration 639 ± 94 ug/uL), as well as good RNA purity and 
integrity (RNA Integrity Number mean 9.1 ± 0.4) (Table 1), 
though there was evidence of some mild contamination 
from purification reagents (260/230 < 2). qPCR was per- 
formed on all RNA preparations to confirm DNA elimin- 
ation; quantification cycle numbers were at least 20 cycles 
lower for RNA samples than equivalent amounts of DNA. 

The modified protocol for RNA and DNA extraction, 
namely the inclusion of lysostaphin (and Dispersin B for 
biofilm samples) in the initial lysis buffer resulted in sub- 
stantial and reproducible increases in RNA and DNA 
yield during harvesting (results not shown). Dispersin B 
catalyzes the hydrolysis of polysaccharide intercellular 
adhesin (PIA), a major constituent of the extracellular 
matrix of Staphylococcus spp. biofilms. Incorporating 
this enzyme in the solution to recover the biofilm from a 
surface, and in the lysis buffer is thought to increase 
recovery by enzymatically degrading the extracellular 
matrix of biofilm, releasing adherent bacteria for recov- 
ery, and exposing them to the lysis solution. The ability 



Table 1 RNA quality and recovery 

Sample ID Concentration 260/280 260/230 Total RNA RIN 





(ng/uL) 






recovery (jig) 




A42 log A 


616.2 


2.15 


1.09 


61.6 


9.6 


A42 log B 


771.4 


2.15 


2.27 


//.I 


n.t. 


A42 log C 


805.5 


2.15 


2.27 


80.5 


9.0 


A54 log A 


823.1 


2.14 


2.20 


82.3 


9.4 


A54 log B 


919.4 


2.13 


1.67 


91.9 


n.t. 


A54 log C 


921.6 


2.14 


2.03 


92.2 


9.4 


A42 stat A 


535.4 


2.10 


2.00 


53.5 


9.2 


A42 stat B 


449.4 


2.02 


1.95 


44.9 


n.t. 


A42 stat C 


493.0 


2.05 


1.66 


49.3 


8.8 


A54 stat A 


403.3 


2.05 


2.05 


40.3 


8.5 


A54 stat B 


630.1 


2.13 


1.89 


63.0 


n.t 


A54 stat C 


397.5 


2.05 


1.89 


39.7 


8.8 


A42 BF A 


732.4 


2.12 


2.17 


36.6 


9.4 


A42 BF B 


738.6 


2.11 


2.18 


36.9 


n.t 


A42 BF C 


887.3 


2.12 


2.19 


44.4 


8.4 


A54 BF A 


606.9 


2.11 


1.97 


30.3 


9.5 


A54 BF B 


337.9 


2.06 


1.74 


16.9 


n.t 


A54 BF C 


421.9 


2.03 


1.99 


21.1 


8.9 



(n.t: not tested, log: logarithmic, stat: stationary, BF: biofilm samples). 



of this simple method to yield an adequate quantity and 
quality of DNA from biofilm-embedded bacteria was an 
important finding and will facilitate future studies of gene 
expression in biofilms. 

Amplification specificity and determination of PCR 
efficiency of reference genes 

Quantification cycle (Cq) for each reaction was plotted 
against the log of DNA concentration, and the slopes of 
the curves were used to calculate the PCR efficiency values, 
which ranged from 1.79 to 1.87, with excellent regression 
coefficients (Table 2). Despite the range in melting points 
of the qPCR products (Table 2), all primers amplified well 
under the conditions listed. Melt curve analysis showed a 
single melt curve for each target gene, and DNA agarose 
gel electrophoresis revealed a single peak for each product. 
Sequencing of the qPCR products matched the sequence 
of the desired target in all cases. 

Comparison of qPCR reaction products 

Melt curve analyses and gel electrophoresis were also 
performed on and compared between samples recovered 
from qPCR reactions performed on cDNA prepared 
using random hexamers, single gene specific primer 
and combination gene specific primer reverse transcript re- 
actions; there were no differences in any of the products 
obtained from those obtained from qPCR or genomic 
DNA. 
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Table 2 Reference and target gene efficiency determination 
and amplicon melting points, r 2 values are regression 
coefficients for the curves 



Candidate 
reference genes 


Slope of 
the curve 


r 2 on the 
slope 


PCR 
efficiency 


Melting 
temperature (°C) 


proC 


-3.78 


0.997 


1.84 


80.0 


gyrB 


-3.85 


0.998 


1.82 


79.5 


rpID 


-3.86 


0.998 


1.82 


75.5 


rho 


-3.68 


0.997 


1.87 


81.0 


rpoA 


-3.89 


0.996 


1.81 


76.5 


ftsZ 


-3.85 


0.998 


1.82 


79.0 


recA 


-3.83 


0.998 


1.82 


81.0 


sodA 


-3.97 


0.996 


1.79 


79.5 


Target gene 










icaA 


-3.92 


0.995 


1.80 


80.5 



Table 4 Pairwise variation in normalization factor for 
combinations of reference genes 

Number of genes used Pairwise variation of normalization factors 



Stability assessment and validation of reference genes 

Relative expression levels (quantification cycle numbers) 
were entered into the Microsoft Excel (Microsoft Canada, 
Mississauga, ON) visual basic application geNorm [12], 
which calculated stability values (M values) for each 
gene (Table 3). GeNorm was also used to calculate 
normalization factors using combinations of the most 
stably expressed genes, and calculated the pairwise 
variation between these factors to identify the optimum 
number of reference genes to use (Table 4). The three 
genes with the lowest individual M values were used in 
combination (rho, recA, gyrB) for subsequent expres- 
sion analysis. 

While stability of expression of reference genes is 
important, reference genes ideally also have a high PCR 
efficiency. The primers we identified had efficiency ran- 
ging from 1.79 to 1.84 (with a value of 2.0 representing 
100% efficiency), and while the final three genes used for 
normalization did not have the highest efficiencies, the 
values are accounted for during normalization of qPCR 
expression data. 

With respect to the identification of reference genes, 
these results are specifically only applicable to the two 

Table 3 GeNorm gene stability (M) values 

geNorm M value 

0.701 
0.604 
0.769 
0.595 
0.663 
0.848 
0.599 
1.150 



Candidate reference genes 

proC 

gyrB 

rp!D 

rho 

rpoA 

ftsZ 

recA 

sodA 



2 vs 3 


0.119 


3 vs4 


0.096 


4vs 5 


0.082 


5 vs 6 


0.112 


6 vs 7 


0.107 


7 vs 8 


0.135 



strains studied under the three growth phases sampled. 
However, the two strains that were studied comprised the 
two main international MRSP clones [3], suggesting that 
these genes will be suitable for broad studies of MRSP gene 
expression. The excellent stability of expression of these 
genes over the broad range of conditions studied also indi- 
cate that these targets will likely function very well as refer- 
ence genes over a generally broad range of conditions. 

icaA expression 

Expression of icaA was significantly higher in the biofilm 
compared to logarithmic and stationary phases (p = 0.0093 
(A42), p < 0.0001 (A54), Figure 1). Individual Tukey's post- 
hoc p-values for individual comparisons were for logarith- 
mic vs biofilm p = 0.015 (A42); p = 0.0001 (A54) and for 
stationary vs biofilm p = 0.015 (A42); p = 0.0002 (A54). 
There was no difference in expression level between loga- 
rithmic and stationary phases for either strain, individual 
Tukey post-hoc p > 0.9999 (A42), p = 0.8258 (A54). 

These findings are unsurprising given the role of this 
gene in formation of polysaccharide intracellular adhe- 
sin, but it is noteworthy that such a profound alteration 
in expression was detectable. Further investigation of ex- 
pression of this gene under other conditions and on a 
variety of surfaces, as well as studies of expression of 
other biofilm associated genes such as MSCRAMMs. 
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Logarithmic Stationary 
Figure 1 Relative expression of icaA. 
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Conclusion 

Proper development of qPCR assays, including reference 
gene assessment, is a critical quality control step in the 
application of this tool. This study has identified a group of 
S. pseudintermedius reference genes, and used those 
reference genes to demonstrate a significant expression 
change in a biofilm associated gene. This information 
provides a vital background for the performance of gene 
expression studies in this increasingly important veterinary 
pathogen. 

Methods 

Bacterial strains and culture conditions 

Two canine S. pseudintermedius strains were chosen, 
representing the two main international MRSP clones, 
sequence type (ST) 68 (strain A42) and ST71 (strain A54). 
These two strains were isolated from clinical infections, 
and were both previously classified as moderate biofilm 
formers using a polystyrene plate assay (unpublished data). 
Isolates were grown in tryptic soy broth supplemented 
with 1% (w/v) dextrose for all growth conditions. Samples 
were isolated as single colonies from streak plates prepared 
from -80°C freezer stocks on Columbia blood agar (Sheep 
blood), and incubated aerobically in a shaker at 37°C. Loga- 
rithmic phase growth was defined in a preliminary study 
(data not presented) as >1 and <6 hours of growth, and 
OD 60 o > 0.5 and <2.0; stationary phase growth was col- 
lected at >12 (12-24) hours, and OD 600 >2.2. For loga- 
rithmic and stationary samples, 10 9 cells were harvested, 
centrifuged (13,000 x g for 30 seconds) to pellet the 
cells, then immediately processed to recover RNA as de- 
scribed below. Biofilm samples were produced by incu- 
bating glass Erlenmeyer flasks unshaken at 37°C for >48 
(48-54) hours, at which time a visible film was present 
on the flask surface. The media was removed, and the 
flasks were washed twice with phosphate buffered saline 
(PBS, pH 7.4). Biofilm was harvested by incubating the 
flask with 10 ug/mL Dispersin B (Kane Biotech, Winnipeg, 
MB) a biofilm degrading enzyme isolated from Aggregati- 
bacter actinomycetemcomitans, in PBS, shaking the flask 
for five minutes at room temperature (the visible film dis- 
persed producing a turbid solution), and then pelleting by 
centrifugation (13,000 x g for 30 seconds) the cells present 
in the solution. The pellet of cells was immediately proc- 
essed to recover RNA as described below. 

Primer design and determination of PCR efficiency 

Eight candidate reference genes (Table 5) were evaluated 
using the validation software geNorm [12]. These genes 
had been examined in S. epidermidis and S. aureus in 
previous studies [13,14], and analogous sequences were 
identified in S. pseudintermedius for evaluation in this 
study. Primers were designed using a combination of 
GeneRunner software version 3.05 (Hasting Software, Inc.) 



and the National Centre for Biotechnology Information 
online primer designing tool (http://www.ncbi.nlm.nih. 
gov/ tools/primer-blast/) using gene sequences for all 
strains of Staphylococcus pseudintermedius (strains ED99 
and HKU10-03) available from Genbank (http://www.ncbi. 
nlm.nih.gov/genbank/). 

Primers for icaA were available from a previous study 
[15] and were validated for use in qPCR as described for 
the potential reference gene primers. 

Pooled chromosomal DNA from the two strains under 
investigation (mixture of equal concentrations from each 
strain) was used to generate dilution series for PCR 
efficiency calculations and for amplicon melting point 
assessment. For efficiency determination, a six fold dilu- 
tion series was prepared in triplicate. Twenty uL reac- 
tions containing from 100 ng to 0.1 pg of total DNA 
were prepared in 96 well plates (Bio-Rad) using the 
LightCycler 480 SYBR Green I Master qPCR reaction 
mixture (Roche Applied Science, Indianapolis, IN), fol- 
lowing the manufacturer's instructions. Primer concen- 
tration was 0.5 uM for each primer in the final reaction. 
A Bio-Rad C1000 thermal cycler with a CFX96 Real- 
Time System and Bio-Rad CFX Manager 2.0 software 
(Bio-Rad Life Sciences, Mississauga, ON) was used to run 
the following optimized thermocycling parameters: de- 
naturation at 95°C for 10 minutes followed by 50 cycles of 
10 seconds at 95°C, 30 seconds at 62°C, 30 seconds at 72°C, 
then a melting curve analysis running from 65 to 95°C with 
steps and measurements every 0.5°C The thermocycler 
software calculated quantification cycle (Cq), efficiency and 
regression coefficients from the recovered data. Amplicon 
identity and primer specificity was confirmed by sequen- 
cing of the PCR products (fluorescent capillary Sanger 
method, Macrogen, Seoul, Korea), melt curve analysis and 
gel electrophoresis. qPCR reactions were tested over a 
range of primer concentrations and annealing temperatures 
and times to determine optimum reaction conditions. 

DNA isolation 

Cells were collected from a logarithmic phase sample 
as described above, and processed using the High Pure 
PCR Template Preparation Kit (Roche Applied Science, 
Indianapolis, IN) following the manufacturer's instruc- 
tions with the exception that lysostaphin (Sigma- Aldrich, 
Oakville, ON) was added at 0.1 mg/mL in the initial lysis 
step. Concentration and purity at 260/280 nm was mea- 
sured using a Nanodrop ND-100 spectrophotometer 
(Nanodrop Technologies Inc., Wilmington, DE), and the 
samples were qualitatively examined for shearing by DNA 
gel electrophoresis. 

RNA isolation 

Each pellet of cells was resuspended in 2 mL PBS (pH 7.4), 
centrifuged to a pellet and the supernatant discarded. The 
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Table 5 Candidate reference and target genes, primers and amplicon sizes 



Candidate reference genes 



Function 



Primer sequence (5'-3') 



Expected size (bp) 



proC 


Pyrrolidine-5-carboxylate reductase 


proC-F gccgaatacaaatgcgcacg 
proC-R aaaaatgcagggccacttcc 


180 


gyrB 


DNA gyrase B subunit 


gyrB-F gcgtccgttgattgaagcg 
gyrB-R aacgtcacttgcaacatcgc 


240 


rpID 


50S ribosomal protein L4 


rplD-F gcctaagaaaatgcgtcg 
rplD-R ccttctggtgttgtgattg 


237 


rho 


Transcription termination factor Rho 


rho-F cacgtaaaagttgctgaattg 
rho-R cctgcttcgatatttctgg 


215 


rpoA 


DNA-directed RNA polymerase sigma factor 


rpoA-F ctatcatcattaccaggtgc 
rpoA-R caaaatttcaacatcactgtc 


231 


ftsZ 


Cell division protein ftsZ 


ffsZ-F gtccattcagtttcgaagg 
feZ-R catgattgttttaacgtcagc 


254 


recA 


Recombinase A 


recA-F gcattaggtgtagatattgataac 
recA-R ggctgcagaaagtttacgc 


228 


sodA 


Superoxide dismutase 


sodA-F cgcaaacttagacagcgtacc 
sodA-R caacaagccaagcccaacc 


227 


Target Gene 








icaA 


N-acetylglucosaminyl-transferase 


icaA-F ttgcccaccttgtgcccacc 
icaA-R tgaggctgtagggcgttggga 


178 



pelleted cells were subsequently processed using the 
Qiagen RNeasy Mini kit (Qiagen, Germantown, MD), 
following the manufacturer's instructions, with the fol- 
lowing changes. The cells were initially resuspended in 
tris-EDTA buffer (pH 7.5) containing 15 mg/ml lysozyme 
(Sigma- Aldrich), 2 mg/mL proteinase K (Sigma-Aldrich), 
25 ug/mL Dispersin B (Kane Biotech Inc., Winnipeg, 
MB) and 0.1 mg/mL lysostaphin (Sigma-Aldrich). All 
optional steps for additional washes or spins during the 
procedure were performed. Logarithmic and stationary 
phase samples were eluted from the purification column 
with two volumes of 50 uL sterile nuclease free water, 
biofilm samples were eluted with one volume of 50 uL 
(volumes chosen so as to obtain final RNA concentra- 
tions above 400 ng/uL). All samples were treated with 
the DNA/ree DNAse kit (Ambion, Austin, TX) as per the 
manufacturer's instructions using a total of 2 units per 
reaction, added in two aliquots. A Nanodrop ND-100 
spectrophotometer (Nanodrop Technologies Inc.) was 
used to measure RNA concentration and purity at 260/ 
280 nm. Three samples of each RNA sample were mea- 
sured and the measured concentrations were averaged. 
Representative samples were submitted for Bioanalyzer 
analysis (Agilent Technologies, Santa Clara, CA). Real- 
time PCR was performed using the efficiency determin- 
ation protocol on RNA samples to confirm the absence 
of DNA using gyrB primers. RNA samples were stored 
at -20 to -80°C until further use. 



Reverse transcription/cDNA preparation 

Two ug of total RNA was used in each 20 uL gene spe- 
cific reverse transcription reaction using the Omniscript 
RT PCR kit (Qiagen), using 10 uM of each primer, 
0.5 mM dNTP, 1 U reverse transcription enzyme, pro- 
vided buffer diluted to lx and the remainder as water. A 
mixture of forward and reverse primers for all nine 
genes was used. Reactions were incubated at 42°C for 
60 minutes. Reverse transcription products were purified 
using the QIAquick PCR purification kit (Qiagen) as per 
the manufacturer's instructions, eluted with 100 uL Qiagen 
EB buffer, and stored at -20°C until further use. 

qPCR for gene stability and icaA 

The same protocol for efficiency testing was used to 
examine stability of potential reference gene and icaA 
expression, with the only change being that only 30 cycles 
were performed. In triplicate, cDNA produced from all 
three growth phases were subjected to qPCR (as described 
for the efficiency evaluation) for each candidate gene. 
In each 20 uL reaction, 5 uL of purified RT reaction 
product (equivalent to 100 ng of RNA prior to reverse 
transcription) was used. No template (negative; water) 
and DNA (positive; 100 ng DNA) controls were included 
in each run. No amplification was identified in negative 
controls. Intra-plate normalization was performed using 
the measured level of icaA in the positive control on each 
plate. 
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icaA expression, statistical analysis 

Expression of icaA was compared between the three 
growth phases for the two strains using the combination 
of three reference genes. Normalization factors, calibra- 
tion factors and relative expression was calculated as per 
Hellemans et al. [16]. A one-way ANOVA with Tukey's 
post-hoc test was perfomed to compare expression be- 
tween growth phases. P-values < 0.05 were considered 
significant. 
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